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ABSTRACT 

The Narrow-line Seyfert I galaxy, 1H0707-495, has been well observed in the 0.3- 
10 keV band, revealing a dramatic drop in flux in the iron Ka band, a strong soft 
excess, and short timescale reverberation lags associated with these spectral features. 
In this paper, we present the first results of a deep 250 ks NuSTAR observation 
of 1H0707-495, which includes the first sensitive observations above 10 keV. Even 
though the NuSTAR observations caught the source in an extreme low-flux state, the 
Compton hump is still significantly detected. NuSTAR , with its high effective area 
above 7 keV, clearly detects the drop in flux in the iron Ka band, and by comparing 
these observations with archival XMM-Newton observations, we find that the energy 
of this drop increases with increasing flux. We discuss possible explanations for this, 
the most likely of which is that the drop in flux is the blue wing of the relativistic.ally 
broadened iron Ka emission line. When the flux is low, the coronal source height is 
low, thus enhancing the most gravitationally redshifted emission. 

Key words: black hole physics - galaxies: active - X-rays: galaxies - galaxy: indi¬ 
vidual: 1H0707-495. 


1 INTRODUCTION 


Sinc e the discovery of t he broad iron K line twenty years 
ago dTanaka et al.lfl995l ). there has been debate as to the 
origin of this distinctive spectral feature. Two physical inter¬ 
pretations fit the time-integrated energy spectrum equally 
well: relativistic reflection, where the iron Ka emission line 
is produced in the inner accretion disc and is asymmetri¬ 
cally broadened due to special and general relativistic ef¬ 
fects in the stro ng potential well of the central supermassive 
black hole (e.g. Fabian et al.ll 19891 : iBrenneman fe Revnoldsl 
l2006l : iDauser et al. 12012k and partial covering absorption, 
where the primary emission is partially covered by dis- 
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tant clumpy clouds or outflowing winds, and this absorp- 
tion pattern mimics the shape of the broad iron line (e. g. 
llnoue fe Matsumotoll2003i ; Miller. Turner, fe ReevesIkOQgk 


This debate comes to a head in the spectral fitting of 
Narrow-line Seyfert I galaxy, 1H0707-495 (z = 0.0411). The 
time-in tegrated spectrum shows a dramatic drop in flux at 
7 keV dBoller et al.ll2002t), which suggests in both the rel¬ 
ativistic reflection mod^jFhbian^etalj 20041 ) and the par¬ 
tial covering model dTanaka et alj 12004 ) that there is an 
overabundance of iron in the emitting/absorbing material 
(unless the absorbing ma terial is outflowing with velocity 
> 0.3c; iDone et al.l 120071 ). In the case of relativistic reflec¬ 
tion, overabundant iron leads to a strong p rediction: the 
broad ened iron L emission line at ~ 0.9 keV. iFabian et aid 
d2009l) showed in a 500 ks XMM-Newton that the residu- 
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als to a phenomenological continuum model show both the 
b roaden ed iron K an d iron L emi s sion l ines. In response, 
iMizumoto. Ebisawa. fe Sameshimal (120141 ') claim that these 
data can be well-fit by a partial covering by a Compton 
thin absorber that envelops a Compton thick absorber, 
however the major features in the spectrum (namely the 
iron-L and iron-K features) are fitted by ad-hoc absorp¬ 
tion edges. They say: “In addition to the cold/thic.k and 
hot/thin absorbers that respectively explain the iron K- 
and L- edges, we had to introduce the ad-hoc edge com- 
ponen t s to explain t he extremely strong edge features.” 
IMizumoto. Ebisawa. fc Sameshimal ( 2014 ) do not discuss 
the implications of adding the edges on the rest of the spec¬ 
tral model, nor do they discuss why the iron-L and iron-K 
edge energies vary independently if they are part of the same 
absorbing structure. The high-r esolution XMM-Newton - 
RGS data were also examined, and lBlustin fc Fabianl J20091 ') 
found no evidence for the narrow emission or absorption 
lines expected from a partial covering modeQ. 

Fortunately, 1H0707-495, in addition to its striking 
spectral features, is also one of the most highly variable 
Seyfert galaxies, commonly known to increase by an order of 
magnitude in flux on timescales as short as one hour. This 
allows for extensive spectral-timing analysis, which can help 
break the degeneracies of traditional spectral modelling. In 
addition to r e vealin g the iron L feature in th e soft excess, 
iFabian et al.l d2009l 'l and IZoghbi et al.1 l|2010l) showed that 
the high-frequency variations in the soft excess from 0.3- 
1 keV lag behind those in the 1-4 keV band by ~ 30 s. These 
authors interpreted this short timescale lag as the reverber¬ 
ation time delay between the primary-emitting corona and 
the inner accretion disc. This short time delay would put 
the corona at a height of < 10r g from the accretion disc 
(using the black hole mass estimate of 2 x 10 6 tfp from th e 
Mb H-X-ray excess-variance relation; iBian fe Zhadl2003lj . 
Later, as predicted by the relativistic reflection model, a 
short timescale lag was also found in the iron K band us¬ 
ing a t otal of 1.3 Ms of archival XMM-Newton data on the 
source dKara et al. Uni). This timing analysis showed that 
the high-frequency variability in the soft excess is correlated 
with that of the iron K band, and that both of these bands 
lag behind the 1-3 keV band by tens of seconds. In other 
words, the energy spectrum of the high-frequency lags re¬ 
sembles the relativistic reflection energy spectrum. To date, 
the simultaneous soft excess and iron K lag has not been 

exp lained with a partia l covering model __ 

iMiller et al.l J20ld l and (and later iLegg et al.ll2012l . for 
a similar source Ark 564) presented a partial covering timing 
model in which longer, low-frequency lags were interpreted 
as reverberation between a compact continuum source and 
distant circumnuclear material. In this model, the amplitude 
of the low-frequency lag increases with energy because the 
continuum photons travel through a medium whose opac¬ 
ity decreases with increasing energy, and therefore the frac¬ 
tion of scattered photons increases with energy. We note 


1 IMizumoto, Ebisawa. fe Sameshimal 1120141 ) also examined the 
same RGS dataset, and claim that there are narrow absorption 
features present. However, these authors do not include a statis¬ 
tical comparison between a featureless continuum and a model 
with absorption lines, so the significance of the putative narrow 
features cannot be accessed. 


that this partial covering model requires the corona to be 
compact because the high-frequency variability must be in¬ 
trinsic to the primary emission. In this model, the high- 
frequency time lags are simply a phase-wrapping artefact 
from the low-frequency lags. If this is the case, then the 
high-frequency lags would be a mirror-image of the low- 
frequency lags (just with a smaller amplitude), however, as 
discussed in the previous paragraph, the high-frequency lags 
sho w a distinct signature in the iron K band (for more d etail, 
see IZoghbi. Uttlev. &; Fabianl l201ll : iKara et al.l 120131 . and 
a complime ntary arg u men t bas ed on lag measurements of 
NGC 6814 in lWalton et al.f2013i b The high-frequency iron K 


lag has now been well-established in a number of source 


(e.g. Zoghbi et al. 

2012|; 

Kara et al. 2013 

; Zoghbi et al. 20lc| 

Alston et al. 2014! 

and 

see Uttlev et al. 

20141 for a review). 


Until partial covering can explain both the high- and low- 
frequency spectral-timing results that have now been widely 
observed, it cannot be considered a viable physical model. 


iGardner fc Pon3 d2014l ) recently presented an alter¬ 
ation to the reverberation lags in 1H0707-495, in which 
reverberation off the accretion flow around a non-spinning 
black hole is diluted by an occulting failed wind at 20 r g . 
This model does require inner disc reflection, but the reflec¬ 
tion fraction is very low because there is no emission from 
within 6r g or more. The authors do not explain the strong, 
very broad iron K lag that is present in this source and many 
others. 


Despite the extraordinary spectral features and vari¬ 
ability patterns of 1H0707-495, it is a relatively faint 
source with a very steep spectrum, which means that it 
has never been detected above 10 keV until now (e.g. 
IWalton et al.l 120131 1. The Nuclear Spectr oscopic Telescope 
Array or NuSTAR dHarrison et al.l l2013l l has the ability 
to make this detection because it carries the first high- 
energy focusing X-ray telescope in orbit, and is 100 times 
more sensitive than previous instruments in the 10-80 keV 
band. Above 10 keV, photons are not photoelectrically 
absorbed and reprocessed. Rather they Compton down- 
scatter, which produces an excess of photons in the 10- 
80 keV range. We call this hard excess the Compton 
‘hump’, and it has been confirmed by NuSTA R in a num¬ 
ber of objects, in cluding NGC 1 365 dR isaliti et aid 120131: 
IWalton et al.ll2014h. MCG-6- 30-15 dMarinucci et al. 2014a 1. 
and Mrk 335 dParker et abl 1 20141 ). NuSTAR even detects 
the associated high-frequency Compton h ump reverberation 
lag in the bright sources, MCG-5-23-16 dZo ghbi et al .lEfflTi] ) 
and SWIFT J2127.4+5654 dKara et alll2015ll . which con¬ 
firms that the variability in the iron K and Compton hump 
bands are correlated, and that they lag behind the contin¬ 
uum by short time delays. This time delay corresponds to 
light travel distances of < 10r g . 


We present the first results of the NuSTAR observation 
of 1H0707-495, which clearly reveals the sharp drop above 
6.4 keV and Compton hump. Given the recent evidence from 
NuSTAR observations and reverberation measurements that 
support the relativistic reflection interpretation, we will fo¬ 
cus on that model for the remainder of the paper. We fit 
the spectrum with a relativistic reflection model in order to 
explore the consequences of including the high-energy emis¬ 
sion and to check that the model works. We do not make 
attempts to fit the spectrum with a partial covering model, 
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Figure 1 . Unfolded spectrum of 1H0707-495 to a powerlaw model with index 0 and normalization 1. The green, red and blue points 
show Obs 1, Obs 2 and Obs 3 of our joint NuSTAR (combined FPMA and FPMB) and Swift observations. The grey points show two 
archival XMM-Newton observations of one of the highest and lowest flux in order to demonstrate the range of fluxes exhibited by this 
source to date. 


Name 

Date 

obsid 

Exposure (ks) 

FPMA/FPMB count rate (counts/s) 

Obs 1 

2014 May 5 

60001102002 

144 

0.008 ± 0.0003/0.007 ± 0.0003 

Obs 2 

2014 June 10 

60001102004 

49 

0.005 ± 0.0004/0.005 ± 0.0004 

Obs 3 

2014 June 28 

60001102006 

47 

0.01 ± 0.0005/0.009 ± 0.0005 


Table 1. Details on the NuSTAR observations of 1H0707-495. Obs 1 and Obs 3 were accompanied by Swift-XRT snapshots. 


as it is not the aim of this paper to determine which model 
is correct based on the NuSTAR spectrum alone. 

The paper is outlined as followed: in Section [2] we de¬ 
scribe the new observations, and we present the spectral 
results in Section [3] We End that the energy of the blue 
wing of the iron line is dependent on flux in this source, and 
discuss possible reasons and implications for this observa¬ 
tion in Section [4] Finally, we summarize our conclusions in 
Section [5] 


2 OBSERVATIONS AND DATA ANALYSIS 

1H0707-495 was observed with the NuSTAR Focal Plane 
Module A (FPMA) and Focal Plane Module B (FPMB) on 
three separate occasions for a total exposure of ~ 250 ks. 
See Table [l] for dates, obsids, exposures and background- 
subtracted count rates. The Level 1 data products were 
processed with the NuSTAR Data Analysis Software (nus- 
tardas v. 1.4.1), and the cleaned Level 2 event files were 
produced and calibrated with the standard filtering criteria 


using the nupipeline task and CALDB version 20140715. 
We chose a circular source region of radius 40 arcsec and a 
background region of radius 85 arcsec. The regions were the 
same for both detectors. While there was spectral variability 
between the three observations, there is little spectral vari¬ 
ability within each observation. The three pairs of spectra 
were binned in order to oversample the instrumental resolu¬ 
tion by at least a factor of 2.5 and to have a signal-to-noise 
ratio of greater than 3<r in each spectral bin. 

Swift/X RT took simultaneous snapshots during Obs 1 
(obsid 00080720001 and 00080720002 for 2 ks and 1.1 ks 
respectively) and Obs 3 (obsid 00080720004, for 1.7 ks). 
There was also a scheduled snapshot during Obs 2; how¬ 
ever, this occurred during an anomaly period for Swift, and 
therefore the data could not be used for spectral analysis. 
The Swift, snapshots were processed using the automated 
processor supplied by the UK Swift Science Data Centre at 
the University of Leicester (lEvans et al.ll2009l 'l. 

We compare these new obse rvations w ith archival 
XMM-Newton observations d Jansen et alj l200ll ; 
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Figure 2. The residuals to a phenomenological continuum model, 
revealing strong features at the iron L, iron K and Compton hump 
bands. The data are from the Swift and NuSTAR observations of 
Obs 1 and Obs 3, which were taken in a similar flux state. 


IStriider et al.l 1200il l. 1H0707-495 has been observed 
for a total of ~ 1.5 Msec with XMM-Newton over 15 
separate observations. The d ata reducti on of these archival 
observations is presented in [Kara et alj (l2013l f. All error 
ranges are quoted at the 90 per cent confidence level for 
each parameter (Ay 2 = 2.71). 


3 RESULTS 

All three observations of 1H0707-495 found the source in 
a low-flux state. Fig. [T] shows the unfolded spectra for our 
three NuSTAR observations along with the corresponding 
Swift snapshots (where possible). Overplotted in grey are 
two archival XMM-Newton observations showing one of the 
highest (obsid: 0511580201) and lowest flux states (obsid: 
0554710801) observed with XMM-Newton. The second of 
the NuSTAR observations was the lowest flux observation 
of the three, and appears to be even lower than the previous 
lowest flux observation taken in 2011 as an XMM-Newton 
Target of Opportunity. Unfortunately, we do not have the 
simultaneous Swift snapshot to confirm that Obs 2 was in a 
similar (or even lower) flux state at soft X-rays. 

NuSTAR, with its high effective area above 6 keV, 
clearly detects a strong drop in the emission around 6-7 keV, 
even when the source is in a low-flux state. Interestingly, the 
energy of this drop is not constant between observations. In 
the high-flux XMM-Newton observation, the drop occurs at 
a higher energy than in the low-flux observations. This vari¬ 
able blue wing was first obse rved in early XMM-Newton ob¬ 
servations llGallo et al.ll2004l f. and we confirm this difference 
now with additional high quality observations. We interpret 
this as a variable blue wing of the iron Ka line, and discuss 
it further in Section [I] 

In Fig. [2] we show the ratio between the combined Obs 
1 and Obs 3 spectrum and the best fit continuum model, 
demonstrating the strong iron L, iron I< and Compton hump 
features. We did not include the Obs 2 NuSTAR spectrum 
because there is no accompanying Swift data. The contin¬ 


uum model was produced by fitting the entire spectrum 
with a black body, powerlaw, two LAOR models at iron K 
and iron L and a Gaussian for the Compton Hump. The 
LAOR lines and Gaussian were then removed, leaving just 
the continuum, consisting of a black body at kT=0.05 keV 
and powerlaw with T = 2.9 with Galactic absorption. Other 
sources observed with NuSTAR have shown similar ratio 
plots with clear iron K and Compton hump features (e.g . 


plots w itn clea r iron l\ and O ompton hump leaturc s (e.g . 
RisahU_jTj^ O|mAjd 2 ] 20 Ul| ^ et al.l I20l4 


Marinucci et al] 2014a 3" Balokovic et al] ' 20141 1. The reflec¬ 

tion features in 1H0707-495 are very pronounced (~ 4 times 
that of the continuum). 

Next, we fit the spectrum with a relativistic reflection 
model in XSPEC. The spectrum can be well described using 
the publicly availa ble relxil l lp, w hich combines the xillver 
reflection grid dCarcfa et al ]H3) with the relconvjp rela¬ 
tivistic convol ution ke rnel for a po int source corona above an 
accretion disc foauser et al.ll2013h . The benefit of the lamp- 
post model is that it returns a physical height to the corona 
and calculates the emissivity profile consistently from the 
source height, rather than assuming a powerlaw or broken 
powerlaw emissivity profile. We emphasize that the lamp- 
post geometry is an extreme simplification of the corona, 
which is likely to be compact, but obviously not a point 
source. We allow for the reflection fraction (ratio of the 
reflected to primary emission, F r /F p ) to vary because the 
self consistent reflection fraction from a point source corona 
will always over predict the true reflection fraction from a 
slightly extended corona. 

We start by fitting the three NuSTAR observations indi¬ 
vidually with the relxilllp model, tying together physical 
parameters, such as the spin, inclination, and iron abun¬ 
dance. The goal was to see which parameters were changing 
between our three observations. Unfortunately, due to the 
quality of the data, we cannot constrain all parameters indi¬ 
vidually. Instead, we start by fitting all three spectra simul¬ 
taneously, tying all parameters together except the normal¬ 
ization. This results in a good fit x 2 /d°f = 251/233 = 1.08. 
We use this as our baseline model, and test untying individ¬ 
ual parameters that could change between our observations 
to see if this improves the fit. The fit improves by letting the 
ionization parameter vary between observations. This results 
in a Ax 2 = 10 for two additional parameters (i.e. the fit is 
preferred with > 99% significance). The ionization param¬ 
eter is constrained to be lower for Obs 1 and Obs 3 than 
in Obs 2: log(£i) < 0.9 erg cm s -1 and log(£s) = 0.4 ± 0.3 
for Obs 1 and Obs 3 versus log(£s) = 1.75lQ 2 g for Obs 2. 
The best overall fit comes from freeing the photon index T 
between the different observations (x 2 = 234/231). The im¬ 
provement is Ax 2 = 17 for two additional parameters. Ti 
and F 3 are consistent within error, and F 2 is significantly 
harder, which is consis tent wit h the observ ed flux and pho¬ 
ton index correlation liWilkins et al]|2014| j. See Fig. [3] and 
Table [2] for the best fit relxilllp model with normalization 
and photon index left free to vary between observations. 
When we let only the source height vary between observa¬ 
tions, the fit was slightly improved from our baseline model, 
but we could not find differences in the source height, as all 
were constrained to be less than 3r g within the 90% confi¬ 
dence limit. 

To confirm that our model was appropriate for the soft 
band as well, we included the simultaneous Swift snapshots, 
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Figure 3. The best fit relativistic reflection model, relxilllp to 
the three NuSTAR observations. The top panel shows the spec¬ 
trum in counts, and the inset shows the unfolded model spectrum. 
The bottom panel shows the x 2 residuals. See Table [2 for more 
details on the fit. 

and added Galactic absorption (IVh = 5.8 x 10 20 cm -2 ) 
and intrinsic absorption to the model. This resulted in a 
X 2 /dof = 288/251 = 1.15. All parameters were the same 
within error, though the ionization parameter was slightly 
higher, log(£) = 2 ± 0.3 erg cm s _1 . The reflection frac¬ 
tion was constrained to be over 15. The higher ionization is 
consistent with previ ous spectral studies of 1 H0707-495 in 
the 0.3-10 keV band. iFabian et aD ll2012l l and lDauser et ahl 
(2012) found that the best fit required two reflectors with 
all parameters the same but with different ionization pa¬ 
rameters. The low-ionization reflection contributes most to 
the iron line, while the high-ionization dominates the soft ex¬ 
cess. The double reflector is an empirical fit to the data, and 
is not motivated by accretion disc models. In these papers, 
the double reflector was interpreted as a patchy accretion 
disc, where the density fluctuates throughout, but it could 
also be interpreted, for example, as an accretion disc that 
has a highly ionized surface layer, with a less ionized layer 
underneath. 

Given the small source height and large reflection frac¬ 
tion found in these observations, we know that the disc ir¬ 
radiation must be very strong. However, the ionization pa¬ 
rameter is constrained to be log(£) ~ 0.3 erg cm s -1 in the 
NuSTAR band and up to log(£)=2 erg cm s~ 4 when includ¬ 
ing the soft excess, both of which are relatively low. The 
ionization parameter £ = 47r Fx/n e , where n e is the electron 
number density in the disc, and Fx is the irradiating X-ray 
flux from T 1000 Ry. If Fx is very high and the ionization pa¬ 
rameter £ is low, this suggests that the density of the disc is 
very high in this object. As an exercise, we compute a rough 
density of the disc given the luminosity and ionization pa¬ 
rameter suggested from our fit to the 0.3-30 keV band. The 


Parameter 

Value 

h (ig) 

2 14 +0,31 
-0.09 

a 

0.97 ±0.03 

in 

43.41 H 

r °ut (r g )“ 

400.0 

IT 

3.05j;°- 25 

r 2 

9 0 7+0.28 

Z,CW -0.17 

r 3 

2 -80io.05 

log(£) (log(erg cm s -1 )) 

0-3±g;l 

A-Fe 

> 7.2 

Ecnt (keV)“ 

300.0 

R 

> 5.8 

AT x 10“ 4 

3.6 ± 1.4 

N 2 x 10 “ 4 

U, ° -0.18 

n 3 x 10- 4 

2 0 +1 ' 6 
z.U—o 2 

X 2 /dof 

234/231 = 1.02 


“ frozen 


Table 2. Best-fit spectral parameters of relxilllp model, cor¬ 
responding to Fig. [3] The differences between Obs 1, Obs 2 and 
Obs 3 can be modelled as differences in just the normalization 
and primary continuum photon index. Therefore N and T are 
denoted with underscores referring to the individual observa¬ 
tions. The normalization is defined as usual in XSPEC, photons 
keV _1 cm _2 s _1 at 1 keV. 

average 1 1000 Ry luminosity from our three observations is 
1.2 x 10 43 erg s _1 . Some of this luminosity is primary emis¬ 
sion, and some is reflection. We are interested in the amount 
of flux irradiating the disc. If we assume that the fraction of 
reflected flux (R > 15) is equal to the fraction of the coronal 
emission that is irradiating the disc, then 15/16 of the total 
luminosity is irradiating the disc. We convert this luminos¬ 
ity to a flux by dividing by 4nr 2 , where r is the distance 
from the corona to the disc, which we find to b e ~ 2r g or 
6 _x 10 11 cm for a 2 x 10® solar mass black hole (jBian fc Zhaol 
l2003h . Given the ionization parameter of 100, we find that 
the electron number density, n e ~ 10 17 cm -3 . The xiller 
and reflionx grids assume reflection of a slab of n e ~ 10 15 
cm~ 3 . The implications of this density discrepancy may be 
important. 

3.1 The low-flux spectrum 

Fig. [l] shows that Obs 2 (red) was found in an extreme low 
flux state, very similar to the one triggered in 2011 with 
XMM-Newton (grey). Unfortunately, there was no Swift 
data available at the time to confirm the soft band similarity 
between these observations. However, we do have simulta¬ 
neous UV coverage with the Swift -UVOT (for the NuSTAR 
observation), and with the Optical Monitor (for the XMM- 
Newton EPIC-pn observation). The U-band magnitudes are 
remarkably similar for these two observations (15.46 ± 0.05 
with Swift and 15.47±0.01 with XMM-Newton), which gives 
us confidence that we are looking at the source in the same 
flux level for the two observations. Furthermore, the rms 
variability of the 3-10 keV band (the common band between 
the two observations) is the same within error. For a 80 ks 
length segment with 3000 s time bins, the XMM-Newton rms 
is 0.48 ± 0.07 and the NuSTAR rms is 0.65 ± 0.28, both of 
which are lower than the rms in liigh-flux observations. The 
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Figure 4. The relativistic reflection model with two different 
ionization parameters fit to the low-flux observations: Obs 2 in 
the NuSTAR band (red), and the 2011 XMM-Newton ToO ob¬ 
servation (grey). The total model is the thick solid line. The 
low- and high-ionization reflectors are the dot-dash and dot¬ 
ted lines, respectively. The model results in a reasonable fit, 
x 2 /dof = 403/355 = 1.14. 


origin of the rapid variability in the low-flux observations 
is likely to include intrinsic variations in the X-ray source 
that cause correlated variations in the reflected emission. 
The similarity of the rms in these two low-flux observations 
lends further support to fitting these two observations si¬ 
multaneously. 

We fit the NuSTAR and XMM-Newton low-state ob¬ 
servations by extending our best fit relxilllp model down 
to the XMM-Newton band from 0.3-10 keV. Again, we in¬ 
clude Galactic and intrinsic absorption into the model, and 
refit. The single absorbed relativistic reflection model to the 
XMM-Newton and NuSTAR low-flux observations gives a 
reasonable fit, x 2 /dof = 470/357 = 1.32, with most of the 
deviations occurring between 1 -2 keV. The statistical fit is 
improved (x 2 /dof = 399/355 = 1.12) if we include a second 
reflector that has all parameters ti ed except the i oniz ation 
parameter and norm alization, as in lFabian et alJ (120121) and 
iDauser et al.l ll 20131) . Details are shown in Fig.[4]and Table[3] 
The high-ionization reflector improves the fit at 1 2 keV. 
The reduced \ 2 is similarly improved if, instead of mod¬ 
elling a patchy disc, we add a separate multi-temperature 
black body component (i.e. diskbb in XSPEC) or a black- 
body with temperature kT = 0.12 keV, both of which fix 
the discrepancy at 1-2 keV. 

With the double reflection model we constrain the low- 
flux state to have a lower source height than the fit to all 
the NuSTAR observations. Interestingly, when we include 
the XMM-Newton low-state spectrum the inclination be¬ 
comes larger, and not well constrained. One possibility for 
this may be because the error bars above 7 keV are large in 
the XMM-Newton spectrum, and increasing the inclination 
helps to broaden the line in this very extreme spectrum. It 


Parameter 

Value 

(Gal) ( cm_2 )“ 

5.8 x 10 20 

Nh ( int ) (cm“ 2 ) 

< 1 x 10 20 


h ( r s ) 


< 1.4 

a 


> 0.988 

i n 


65.0 ±14.0 

r out (fg)“ 


400.0 

r 


2.57 ±0.06 

l°g(£l) (log(erg cm 

s- 1 )) 

3.2 ±0.3 

log(? 2 ) (log(erg cm 

s- 1 )) 

■i q+ 0.03 
± * z -0.1 

^Fe 


> 9.5 

Ecut (keV)“ 


300.0 

R 


» 10 

Ni x10“ 7 


0.1 ±0.4 

N 2 x10- 7 


10. ±0.4 

X 2 /dof 


403/355 = 1.14 


“ frozen 


Table 3. Best-fit spectral parameters of Galactic and intrinsically 
absorbed reflection model (with two different ionization param¬ 
eters) to the low flux XMM-Newton and NuSTAR observations, 
corresponding to Fig. [T] 

may also be that the inclination is increasing because the 
model is trying to fit the excess at 2-3 keV with the broad 
iron line. It is possible that this 2-3 keV excess is not iron, 
but actually overabundant sulphur at ~ 2.3 keV. Unfortu¬ 
nately, we cannot account for overabundant sulphur with 
the current relxilllp model. The powerlaw flux is so low 
in this observation that we cannot constrain it with the cur¬ 
rent data. Therefore we cannot put strict constraints on the 
reflection fraction, though it is clearly very high. 

The self-consistent reflection fraction for a point source 
of h = 1.4 is R = 8.2. However, our source height is an up¬ 
per limit, and therefore that self-consistent reflection frac¬ 
tion is a lower limit, which greatly increases as the source 
height decreases down to the event horizon (e.g. our true 
minimum in x 2 " s P a ce is h = 1.17, which corresponds to a 
reflection fraction of R = 22.3). If, instead of fitting the 
reflection fraction separately, we fit the reflection spectrum 
self-consistently for a point-source above the disc (fixRe- 
HFrac=l in the relxilllp XSPEC model), we obtain a 
fit with y 2 /dof = 536/358 = 1.49. In this scenario, the 
higli-ionization reflector is no longer required by the data, 
however this is still a statistically worse fit that the single 
reflector without a fixed reflection fraction (Ax ' 2 = 63 for 
one less degree of freedom). The fit is much improved with 
an additional blackbody of temperature, kT = 0.12 keV 
(X 2 /dof = 404/356 = 1.14). This tells us that there is some 
discrepancy between the data and the simple model of a 
static point-source above an infinite thin slab of density 
n e = 10 15 cm' 3 . These high-quality data are pushing the 
limits of our theoretical models, and therefore more work 
on the models is required, but is beyond the scope of this 
paper. 


4 DISCUSSION 

We have presented the analysis of the first NuSTAR ob¬ 
servations of 1H0707-495. The source went into an extreme 
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Flux (erg/s/cm 2 ) 



Figure 5. Left.: The rest-frame energy of the blue wing of the iron line versus the 3—10 keV flux. The blue wing energy increases with 
flux up to ~ 6 X 10^ 13 erg/s/cm 2 , above which the line energy scatters around 7.2 keV. Right: The powerlaw photon index from the 
‘continuum’ band (1.5—2.5 keV) versus the flux at 3—10 keV. We do not include our newest observations in this as the continuum cannot 
be well constrained from the short Swift snapshots. We see a similar dependence on flux for the edge energy and the photon index, 
perhaps indicative of an underlying change in the corona with flux. 


low-flux state during these observation, but can still be ob¬ 
served up to 30 keV thanks to the unprecedented effective 
area of NuSTAR in the hard X-ray band. The data are well 
described by a relativistic reflection model, where the lamp- 
post source height is constrained to be ~ 2r g . As in previou s 
spectral studies (e.g. lFabian et al1l2009l ; lDauser et alJl20l3) . 
the spin and iron abundance are very high, a = 0.97 ± 0.03 
and ^Fe > 7.2, respectively. The reflection fraction is also 
ve ry high in these obs ervations, similar to what was found 
in iFabian et al.l d2012h in a triggered XMM-Newton obser¬ 
vation of 1H0707-495 in an extreme low state. The high re¬ 
flection fraction in low-flux observations can be explained 
by the strong gravitat ional light bending effect s from a low 
coronal source height llMiniutti fe Fabianll2004l 'l. A high re- 
flection fraction is also evidence for a high black hole spin 
dDauser et al.ll20l3 : iParker et al .1120141 ). 

In this section, we discuss the apparent variable blue 
wing of the iron line, as seen by eye in Fig. |T] We quantify 
this effect using all 15 archival XMM-Newton observations 
of 1H0707-495 in addition to these recent NuSTAR obser- 
vations. The photon index is observed to change with flux 
dWilkins et al.l l2014l. and therefore it is important to get a 
good description of how the underlying continuum changes. 
To this end, we fit the 1.5-2.5 keV band with a power law for 
each observation. We freeze the photon index and normal¬ 
ization of the underlying continuum, and then examine the 
3-10 keV band (the iron line band). We fit the excess above 
the continuum with an additional powerlaw with an edge. In 
effect, we are fitting the red wing of the iron line as a power 
law, and the blue wing of the line as an edge. This is a good 
diagnostic for measuring the energy of the blue wing. All 
observations were well fit by this model. One 40 ks observa¬ 
tions taken in 2007 (obsid: 0506200301) had too few counts 
to put meaningful constraints on the edge of the iron line, 
so for this observation, we fit it together with the follow¬ 
ing observation (obsid: 0506200201), which was contiguous 


with the first. In Fig. [5]-le/f we plot the best fit edge energy 
and the 3-10 keV flux for the 14 independent XMM-Newton 
observations in black, and our three new NuSTAR observa¬ 
tions in green, red and blue. We find that, especially for flux 
levels below 8 x 1CP 13 erg cm' 2 s -1 , the energy of the blue 
wing increases with increasing flux. The edge energies of the 
two highest flux observations do not follow the increasing 
trend. A similar trend with flux is found in the continuum 
(1.5-2.5) keV powe rlaw photon in dex (Fig. fright). Similar 
to recent results by I Wilkins et all d2014ll . the photon index 
increases with flux. However, again, the trend does not hold 
for the two highest flux observations. 

I Gallo et al.l d2004l l and lFabian et all d2004h wrote back- 
to-back papers on the variable blue wing in 1H0707-495. 
These papers were based on the first two XMM-Newton 
observations of the source. They found that the spec¬ 
tra could be equally well explaine d by a changing par¬ 
tial covering absorber dGallo et al.l [20041 1. or by changes 
in the emissivity profile of the accretion disc with flux 
dFabian et al.ll2003 l. If the data are described by the par¬ 
tial covering absorber, then the change in edge energy is 
best explained by a high speed outflow in the observa¬ 
tion with the higher edge energy. The absorber with vari¬ 
able edge energy from a Compton thick outflow was also 
used recently to describe the X M M-Ne wton observations 
dMizumoto, Ebisawa. fc Sameshimall2014l l. This interpreta¬ 
tion would imply that the flux is increasing as the out¬ 
flow velocity increases, which is possible if we assume that 
no new absorbers are generated and the covering fraction 
decreases with higher velocity outflow. However, clumpy 
Compton thick absorption along our line of sight will result 
in a strong scattered component, unle ss the absorbers satisf y 
very restrictive geometric constraints dRevnoIds et al.ll2009i l. 
No clear narrow features are observed in the XMM-Newton 
EPIC-pn or RGS data, though, of course, this non-detection 
does not rule out the possibility of the specific geometry of 
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multiple absorbers along our line that reabsorb all the scat¬ 
tered emission. 

The other interpretation explored to explain the shift¬ 
ing blue wing was changes in the coronal height in the rel¬ 
ativistic reflection model. If in the low-flux state, the X-ray 
emitting source is compact and the source height is small, 
then only a small portion of the accretion disc will be irra¬ 
diated due to strong gravitational light bending effects, and 
the observer will see mostly the gravitationally redshifted 
emission. The flux increase can be explained as an increase 
in the coronal source height, which then causes larger radii 
to be irradiated. This produces emission that is not grav¬ 
itationally redshifted. In addition to light bending effects, 
the structure of the disc itself could play a role in enhancing 
the gravitationally-redshifted emission at low flux levels. In a 
standard Shakura & Sunyaev accretion disc near its Edding¬ 
ton limit, the scale height of the disc far from the innermost 
stable circular orbit can be as large as 5r g . If the coronal 
source height is less than 5 r g , then, at medium inclina¬ 
tions, it would be eclipsed by the large scale height disc at 
larger radii. The eclipse of the primary emission would cause 
the reflection fraction to increase, as we would still observe 
reflected emission from the far side of the disc (which is not 
eclipsed). In this model, when the source height is higher, 
it will not be eclipsed and so the observed flux increases. 
The large source height means that larger portions of the 
disc are irradiated, and the line centroid is clearly detected. 
Both light bending and the disc eclipse would have the ef¬ 
fect of steepening the emissivity profile with decreasing flux. 
Similar affects have been found for the evoluti on of the emis¬ 
sivity break radius with flux in 1H0707-495 dWilkins et al.l 
I2014h . and the steepening of the inner emissivity index with 
increasing reflecti on fraction in the ‘twin’ to 1H0707-495, 
IRAS 13224-3809 rtChia.ng et aHl20!4li . 

The two highest flux observations do not agree with 
the trend of increasing edge energy with flux. One possi¬ 
ble explanation for this is that the corona is more vertically 
extended in the high flux observations than in the ones be¬ 
low 8 x 10 -13 erg cm -2 s _1 . While more photons can es¬ 
cape the gravitational pull of the central black hole, there 
is still a significant portion of the continuum that is irra¬ 
diating the innermost regions of the disc. The fact that we 
also see this turnover at high fluxes for the photon index is 
al so expected i f the corona is more extended. For example, 
in iKara et aD (l2013l 'l. we showed that the energy spectrum 
of the high-frequency variability is harder than that of the 
low-frequency variability, suggesting that the corona emits 
a harder cont inuum at sm aller radii. X-ray microlensing in 
quasars from IChen et al.1 (1201 ll ) also show that the harder 
X-ray emission originates from smaller radii than the soft 
emission. The spectral hardening of the two highest flux 
observations would then suggest that some fraction of the 
coronal emission in these two observations is produced at 
small radii. We have not fit the archival XMM-Newton ob¬ 
servations with an extended c orona, and therefore we ca nnot 
quantify the extent. However IWilkins fc Fabianl (l2013l l has 
suggested through ray-tracing simulations of an extended 
corona, that 1H0707-495 is extended out to 20-30 r g , though 
most of the reflected emission comes from within 5 r g . 

We also explore the effects of ionization on the broad 
line, but find that this is likely not the dominant effect. As 
the ionization increases, there are Ka transitions from more 


Fe ions, and we see emission from the H- and He- like ions (Fe 
xxv-xxvi) up to ~ 6.9 keV. This effect is only appreciable 
when the ionization parameter is high (above log(£) ~ 3). 
Our spectral fitting shows that the ionization increases with 
decreasing flux (consistent with stronger light bending in 
the low-flux state). If ionization changes were driving the 
change in edge energy, then we would expect the opposite 
correlation with flux. 

Lastly, changes in disc inclination would also change 
the blue wing of the iron line, but this is not likely to be the 
cause. If the disc were precessing, then the edge of the wing 
would decrease with increasing flux. In the precessing disc 
model, as the inclination of the disc becomes smaller(more 
face-on), we see a higher flux because of the solid angle in¬ 
crease. However, this would also decrease the Doppler boost¬ 
ing effects, and so the line would appear to be less broad (i.e. 
the energy of the blue wing of the line would decrease). Fur¬ 
thermore, the second and third NuSTAR observations were 
taken days apart from each other, and the disc is unlikely 
to change on such short timescales. 

The variable blue wing energy complicates our spectral 
modelling. If the changing emissivity profile affects the en¬ 
ergy of the blue wing of the line, then we cannot use the 
blue wing as a clear indicator of the inclination of the accre¬ 
tion disc. Usually, the blue wing energy is constant because 
it is determined by how much Doppler shifting is occurring. 
Therefore, it is likely that by fitting the NuSTAR data alone, 
we are inferring a lower disc inclination because it appears 
that there is less Doppler broadening. 

In order to account for this variable blue wing in our 
spectral modelling, we fit the NuSTAR observations with the 
same relxilllp model, but this time freeze t he inclination 
to 55.7 degrees. This is the inclination found in lFabian et al.l 
(l2009ll from fitting a relativistic reflection model to 500 ks of 
XMM-Newton data when the 3-10 keV flux was ~6x 10 -13 
erg cm -2 s _1 . For this test, we take 55.7 degrees to be the 
‘true’ inclination, and use it to fit the recent NuSTAR spec¬ 
tra. This model gives a good fit, \ 2 /dof = 265/230 = 1.15, 
and other parameters change slightly to account for this 
> 10 degree difference in inclination. The effect of increas¬ 
ing the inclination is generally to decrease the height of the 
source. This, in return, increases the ionization and possibly 
the reflection fraction, which is to be expected as stronger 
gravitational light bending occurs when the source is lower. 
When the inclination is fixed at 55.7 degrees, the height is 
more tightly constrained to be even lower, h < 1.5 r s . The 
ionization increases to log(£) = 1.6 ± 0.7 erg cm s , and 
the reflection fraction is constrained to be slightly higher, 
R > 7.8. All other parameters remain largely the same. The 
spin is still high, a = 0.96 ± 0.01, and the iron abundance is 
still high AFe > 8. 

We can track the changes in the energy of the blue wing 
of the iron line in 1H0707-495 because it has been observed 
for several years in many different flux states, but in addition 
to the impressive coverage, 1H0707-495 is special because it 
has a high iron abundance, which makes it easier to measure 
the iron features. The question of the high iron abundance 
in 1H0707-495 has long been a criticism of the reflection 
modelling in this source. Our new observations also require 
a high iron abundance. In fact, if we ignore the NuSTAR 
data below 7.5 keV (i.e. below the iron line) and fit the high 
energy spectrum, we find that the fit requires an iron abun- 
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dance greater than 2.9 times solar abundance at the 68% 
significance. While this is not extremely over-abundant or 
highly significant, it does suggest that the high iron abun¬ 
dance is not motivated by the iron K edge alone, but also 
by the amount of photoelectric absorption in the 8-13 keV 
band. 


5 CONCLUSIONS 

We have presented the spectral analysis of the first NuSTAR 
observations 1H0707-495. Our main findings are: 

• NuSTAR confirms a strong drop in flux at the iron K 
band, and observes the Compton hump in this object for the 
first time. 

• The NuSTAR observations were taken when the source 
was at a relatively low-flux. The lowest flux observation 
(Obs 2), appears to be at an extreme low-flux state. 

• The three NuSTAR observations can be well fitted by 
a relxilllp model with high spin and a high iron abun¬ 
dance. The iron abundance derived from just fitting above 
7.5 keV also suggests a super-solar iron abundance with 68% 
significance. 

• The lowest flux observation, Obs 2, can be well fit 
with the extreme low-flux XMM-Newton observation with 
a model of relativistic reflection off a patchy (non-uniform 
density) disc. 

• The energy of the blue wing appears to change with flux 
in this source. The mostly likely explanation for this depen¬ 
dence on flux is that at low fluxes, the coronal source height 
is lower, and therefore we observe the more gravitationally 
redshifted emission at lower flux. If we assume that the ‘true’ 
inclination is the one found when the source was in a high- 
flux state and fit our NuSTAR spectrum again fixing the 
inclination to a higher value, then the source height for our 
new observations is constrained to be less than h < 1.5r g . 
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